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v  20  ABSTRACT  (ConHnumd) 

1  vides  plasma  diagnostics  capabilities  by  calculating  emission  intensities  of  the  second 
positive  and  the  first  negative  band  systems  of  N£  and  NSj*. 

The  code  is  applied  to  breakdown  calculations  in  N2  and  the  results  are  in  reasonable 
accord  with  experimental  data.  The  wave  optics  code  for  energy  deposition  is  also 
coupled  to  an  air  chemistry  code  (CHMAIR)  to  describe  the  breakdown  off  a  reflecting 
surface.  ~j 

An  insight  is  provided  into  the  plasma  absorption  by  using  a  simple  slab  density  pro¬ 
file  for  the  electrons.  This  models  roughly  the  profile  expected  in  air  breakdown  off  a 
reflecting  surface.  Simulation  results  show  a  steep  localized  electron  density  profile  at 
the  interference  maximum.  ^These  simple  calculations  show  the  plasma  reflectivity  to  be 
small  (<  25%)  for  1  cm  lighted  high  electron  densities,  Ne  1015  cm*3  >  Nc  =  1013 
cm*3.  Simulation  results  also  show  the  possibility  of  attaining  high  electron  density. 
However,  in  breakdown  off  a  reflecting  surface,  the  microwave  is  decoupled  from  the 
plasma  when  the  absorption  length  is  less  than  one-half  the  wavelength  of  the  light. 
Important  questions  to  be  addressed  in  the  future  include  the  possibility  of  plasma 
maintenance,  gas  heating,  shock  formation,  plasma  decoupling  and  effects  of  strong  non¬ 
linear  forces. 
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I .  INTRODUCTION 

The  interactions  of  the  high  power  microwave  radiation  with  air, 
matter  and  plasmas  have  potential  for  various  applications,  and  hence  merit 
detailed  studies  and  analysis.  A  detailed  study,  theoretical  and 
experimental,  is  necessary  in  order  to  delineate  possible  regimes  for 
applications  purposes. 

Studies  can  be  performed  analytically  to  provide  answers  in  the 
limiting  cases.  However,  the  phenomena  of  the  microwave  interaction  with 
air,  matter  and  plasmas  are  complex  and  require  time  and  space  dependent 
treatments  of  a  large  number  of  physical  processes.  This  can  be  carried 
out  in  a  self  consistent  manner  by  developing  detailed  computer  codes  which 
treat  part  or  all  the  phenomena. 

This  report  deals  with  and  describes  computer  codes,  developed  at  NRL, 
for  the  analysis  of  the  pulsed  microwave  interaction  with  nitrogen  (N2) . 

The  interaction  is  considered  from  the  first  principles  with  the  microwave 
absorption  by  one  or  few  free  electrons,  present  in  the  gas,  the  electron 
build  up,  the  breakdown,  and  the  microwave  absorption  and  reflection  by  the 
plasma  electrons.  The  computer  code  which  describes  this  phenomena  is 
basically,  consists  of  two  codes  coupled  together.  One,  describes  the 
electromagnetic  radiation  (E  and  M  Code)  and  the  other  describes  the  physics 
and  chemistry  of  the  ionization  and  the  deionization  in  N2  (CHEM-N2)  Code. 
The  E  and  M  Code,  however,  can  also  be  coupled  into  an  air  chemistry  code  to 
describe  the  breakdown  in  air,  as  discussed  in  the  report. 
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In  Section  1  the  E  and  M  Code  and  the  solution  of  the  wave  equation  is 
described.  Section  2  gives  the  description  of  the  pure  N2  chemistry  code, 
and  Section  3  deals  with  the  breakdown  calculations  in  N2,  where  the  break¬ 
down  power  thresholds  are  compared  with  experimental  data  for  several  nitro¬ 
gen  pressures.  In  Section  the  E  and  M  Code  coupled  to  an  air  chemistry 
code  is  utilized  to  calculate  breakdown  in  air  near  a  surface.  For  air 
breakdown  near  a  surface  one  can  utilize  microwave  powers  below  the  threshold 
for  air  breakdown  by  a  factor  of  k  and  still  be  able  to  break  down  the  air. 
This  is  due  to  the  fact  that  the  microwave  reflection  from  the  surface 
interferes  constructively  with  the  outgoing  wave  resulting  in  doubling  the 
electric  field  at  a  distance  of  — from  the  surface. 

2.  The  E  and  M  CODE 

In  order  to  describe  the  microwave  absorption  by  a  plasma  one  must  know 
the  electric  field  at  all  positions  in  the  plasma.  To  do  this  we  consider 
the  solution  of  one  dimensional  wave  equation1 
32E  (Ua  3E  32E  CD2 
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for  a  plane  wave  incident  normal  to  a  plasma.  The  electron  neutral  (electron- 
ion)  col liaion  frequency  is  denoted  by  v  and  the  plasma  frequency  is  defined 

C 

as  ,„2» - where  the  notation  is  standard.  Since  the  absorption 

p  m 

length  at  microwave  frequency  is  long  at  the  critical  density  =  m  u)2/Wre2 

we  can  specify  an  essentially  free-space  boundary  condition2  for  the 

incoming  light  as  3E  3E  3 
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where  U^n  is  the  amplitude  of  the  known  incident  light.  This  boundary 
condition  is  also  applicable  to  a  plasma-vacuum  interface  as  constructed 
in  experiments  using  a  mylar  divider.  The  right  side  boundary  condition  in 
the  plasmas  specified  as 
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where  j  «  <rEz  *  Ez*  f°r  l°w  densities,  the  right 

hand  side  of  equation  (2)  is  similar  to  a  free-space  propagating  to  the  right. 
Equation  (1)  is  differenced  implicitly  to  retain  stability  for  large  times 
steps  in  a  high  density  plasma,  m2  2  vc2  +  a)2.  The  implicit  difference  form 
of  equation  (1)  is 

_n+l 


'j 


2E*  +  E1)-1 

j  J 


,n+l 


.n-1 


U)_ 


(At)5 


02  +  V2 


At 


s-.rH-l 


(62E) 


+  2(5aE)  ”  +  (5^) 


n-1 


0) 


u>_ 


—  =  0 


(Ay)* 


u,2  +  v2 


(3) 


where  (5  2E)°  =  Ej+^  -2E”  +  E^  This  can  also  be  cast  into  a 
tri-diagonal  system  of  the  form 
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A  solution  of  a  tri-diagonal 


system  of  equations'1 
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The  solution  technique  is  then  to  specify  the  boundary  condition  and  sweep 
the  mesh  space  for  the  Gj's  ^  Fj's*  A^ter  which  we  sweep  again  to  obtain 
the  updated  electric  field. 

The  stability  analysis  of  the  implicit  scheme  can  be  made  simple  by 
dropping  the  last  term  of  Equation  (1).  This  term  is  typically  not 
important  in  determining  the  stability  criterion.  We  rewrite  Equation  (1) 
without  this  term  as  a  system  of  two  coupled  equations. 
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We  find  after  some  algebra  that 
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where  y  s  cAtsin  (kAv/2)  .  Stability  requires  the  eigenvalues  of  the 
Ay 

amplification  matrix  A  to  be  less  than  or  equal  to  1.  One  can  show  that  the 
eigenvalues  are  bounded  by 
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and  hence  the  implicit  scheme  is  unconditionally  stable  for  arbitrary  or. 

The  boundary  conditions  are  explicitly  differenced  and  are  coupled  to 
an  explicit  form  of  Equation  (1)  to  obtain  boundary  value.  This  explicit 
form  of  Equation  (1)  is: 
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Once  the  boundary  values  are  obtained,  we  use  the  implicit  difference 
equation  for  all  interior  points.  The  boundary  values  are  used  as  "first" 
values  in  sweeping  the  tri-diagonal  system  of  equation.  Since  locally  the 
explicit  differencing  is  stable  for 
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the  boundary  must  be  at  low  density  to  maintain  stability,  a>  2  £  v*  +  m2. 
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In  practice,  the  implicit  scheme  must  satisfy  the  explicit  stability  criterion 
to  maintain  accuracy  and  minimum  phase  error.  In  considering  microwave  break¬ 
down  off  a  reflecting  surface,  this  is  important  since  the  interference 
maximum  determines  the  breakdown  field. 

As  a  check  on  the  code,  we  perform  simulations,  on  free-space,  propaga¬ 
tion  and  reflection  off  metal  surface  in  a  vacuum.  Figure  (la)  shows  the 
reflected  light  for  a  propagating  electromagnetic  wave  in  vacuum  for  a  single 
pass  to  the  reflecting  surface  and  back  to  the  point  of  entry  of  the  radia¬ 
tion.  The  system  length  is  kL  =  15  (where  k  =  2rr/\)  so  the  return  time  of  the 
radiation  is  t  uu  *  30  for  cAt  =  0.9Ay.  The  code  tracks  the  propagation 
velocity  well,  however,  in  the  transient  time  tut  £  30,  the  numerics  has  intro¬ 
duced  errors  in  the  reflected  wave.  The  average  reflected  wave  should  be 
zero  during  the  transient  period,  yet  we  see  the  average  reflected  wave  of 
roughly  10$  of  the  incident  light.  Figure  (lb)  shows  the  spatial  electric 
field  (in  normalized  variables) .  Since  the  launched  incoming  wave  is  unity, 
the  combined  wave  once  reflected  from  the  surface  should  be  of  amplitude  two. 
We  see  the  implicit  scheme  has  introduced  large  errors  in  the  phase  of  the 
wave  and  therefore  the  reconstructed  wave.  Figure  (2a)  shows  a  similar  test 
with  a  reduced  time  step  cAt  *  0.2Ay.  However,  in  this  test,  the  reflected 
wave  in  the  transient  period  is  less  than  1$  of  the  incident  light.  The 
smaller  phase  error  reconstructs  the  interference  patterns  to  its  proper 
maximum  of  two  as  can  be  seen  in  Figure  (2b) .  An  extensive  parameter  study 
shows  good  results  are  obtained  using  the  implicit  scheme  for  cAt  £  0.2Ay. 

Also  as  a  check,  we  consider  light  normally  incident  on  a  linear  density 

profile.  The  density  profile  in  this  collisionless  plasma  (vfi  *  O.luj)  rises 

from  zero  to  the  cut  off  density  N  with  a  density  scale  length  of  kL  =  10. 

c  s 

The  Airy  function  solution  requires  standing  wave  pattern  toward  the  incident 
light  and  evanescent  wave  in  the  overdense  plasma.  Figure  (3)  shows  the 
solution  obtained  for  the  electric  field  using  the  implicit  scheme  with  the 
boundary  conditions.  Analytic  solution  shows  the  computed  field  is  accurate 
in  location  of  the  first  large  maximum,  the  amplitude  of  the  reflected  light 
and  the  amplitude  of  the  largest  maximum.  The  various  test  problems  thus 
confirm  the  validity  of  the  code. 

3.  N2  BREAKDOWN  AND  OPTICAL  DIAGNOSTICS 

To  describe  the  nitrogen  breakdown  by  the  microwave  radiation,  we  have 


vwvx; 


J _ L 


V  V 

1 


20  SO  40  SO 


(a)  (b) 

Fig.  1  —  Test  run  with  cAt  =  0.9Ay  for  a  plane  wave  reflected  from  a  surface  (a)  reflected 
wave  (VOS/VEO)  versus  time  (toj)  and  (b)  electric  field  (VOS/VEO)  versus  position  (ky) 
where  VOS  and  VEO  indicate  the  oscillatory  and  the  thermal  velocities  of  the  electron, 
respectively. 


(a)  (b) 

Fig.  2  —  Test  run  with  cAt  =■  0.2Ay  for  a  plane  wave  reflected  from  a  surface  (a)  reflected 
wave  (VOS/VEO)  versus  time  (tw)  and  (b)  electric  field  (VOS/VEO)  versus  position  (ky) 
where  VOS  and  VEO  indicate  the  oscillatory  and  the  thermal  velocities  of  the  electron, 
respectively. 
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(a)  (b) 

Fig.  3  —  Test  run  tor  a  plane  wave  normally  incident  on  a  linear  plasma  profile. 
The  graph  shows  the  electric  field  ( VOS/VEO)  versus  position  (ky). 
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developed  an  N2  chemistry  code.  The  code  calculates  the  time  developments 

of  N2+(x),  N2+(B),  N4+,  N,  N(2D),  NaCA3!),  N2(B3tt)  ,  N2(C3tt)  ,  N  ,  T  ,  T  ,  and 

e  e  v 

T  ,  where  the  last  three  symbols  designate  the  electron,  the  N2  -  vibrational 
and  the  gas  temperatures,  respectively.  The  purpose  for  the  calculation  of 
the  population  densities  of  several  excited  states  is  to  provide  volume 
emission  rates  for  optical  diagnostics  of  the  plasma  and  the  investigation 
of  the  ionization  from  the  excited  states. 

3.1  N2  CHEMISTRY 

A  free  electron  gains  energy,  in  the  field  of  the  microwave  radiation, 
it  undergoes  both  elastic  and  inelastic  collisions  with  N2,  thereby  losing 
part  or  all  of  the  energy  gained  from  the  field.  When  the  electron  energy 
reaches  and  overshoots  the  ionization  threshold  for  N2  an  ionization  event 
may  occur  according  to  reaction  (11). 

e  +  N2  -»  N2+  +  2e  (11) 

Higher  energy  electrons  could  ionize  N  through  the  dissociative  ionization. 
However,  the  threshold  for  this  process  is  higher  by  ~  10  eV  compared  to  that 
for  reaction  (11).  Therefore,  process  (11)  is  the  predominant  reaction  for 
the  breakdown  ionization  by  microwave  radiation. 

The  deionization  of  the  plasma  occurs  primarily  through  the  dissociative 
recombination  of  the  nitrogen  ion  N2+  and  its  cluster  ion  N4+,  i.e., 

N2+  +  e  -*  N  +  N 

N4  +  e  -4  N2  +  N2 

The  N4+  cluster  forms  via  a  three-body  process4'5  i.e., 

n2+  +  n2  +  n2  -»  N4+  +  N2  (llo 

and  is  obviously  important  at  atmospheric  pressures  and  lower  gas  tempera¬ 
tures.  When  the  gas  temperature  rises,  however,  the  formation  of  N4+  is 
inhibited  by  the  slowing  down  of  the  forward  reaction  in  (14)  and  by  the 
collisional  breakup  of  N4+  (i.e.,  the  reverse  reaction  of  (lV)  becomes 
operative . 


(12) 

(13) 


The  deionization  also  proceeds  via  the  three-body  collisional 
recombination,  however,  it  may  become  important  after  breakdown  and  when 
the  electron  density  is  quite  high.  Finally,  the  electron  loss  occurs 
also  through  diffusion. 

3.2  ELECTRON  ENERGY  GAIN  AND  LOSS  PROCESSES  IN  N2 

The  electron  energy  loss  processes  in  N2  proceed  by  the  elastic  and 

inelastic  collisions.  The  energy  loss  by  elastic  collisions  is  dependent 

on  v  ,  where  v  is  the  momentum  transfer  collision  frequency.  The  collision 
c  c 

frequency  depends  on  the  electron  energy  (temperature)  and  the  momentum 
transfer  cross  section.  The  cross  section6  is  shown  as  a  function  of  the 
electron  energy  in  Figure  (k) .  By  folding  this  cross  section  with  the 
appropriate  electron  velocity  distribution,  one  obtains  the  momentum  transfer 
collision  frequency.  Thus  the  rate  of  the  energy  loss  by  an  electron  through 
the  elastic  collisions  can  be  expressed  as 

2m  3  3 
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where  m  is  the  electron  mass  and  M  is  the  mass  of  the  nitrogen  molecule. 

The  electron  inelastic  energy  loss  processes  in  N2  are  the  excitation, 


the  dissociation  and  the  ionization 

processes  are  expressed  as 
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the  nitrogen  molecule. 
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in  addition  to  Equation  (1).  Equation  (16)  Indicates  the  excitations  of  the 
Ns  ground  state  vibrational  levels,  while  Equation  (17)  represents  the 
excitations  of  the  electronic  states  of  N2.  The  N2  electronic  states  include 
the  triplets  and  the  singlets  whose  cross  sections  are  presented  and  discussed 
elsewhere.7  The  cross  section  for  the  N2  disassociation,  however,  is  a 
composite  of  a  large  number  of  singlet  states  which  predissociate.  Thus,  the 
energy  loss  by  electrons  in  N2  can  be  summed  up  to  include  the  losses  due  to 
the  excitations  of  the  ground  state  vibrational  levels,  the  excitation  of 


the  electronic  states  Aae,  B3e,  C3tt,  a^rr,  the  dissociation,  and  the 
ionization  of  N2» 

Briefly,  the  vibrational  cross  sections  are  from  the  measurements  of 
Schulz8  and  are  well  established  as  to  the  peak  of  the  total  and  the 
shape  of  the  individual  cross  sections8  10 .  The  cross  section  for  the 
ionization  is  from  the  measurement  of  Rapp  and  Golden11  and  the  cross 
section  for  the  dissociation  is  from  the  measurements  of  Winters12  and 
Zipf  and  Mclaughlin13.  The  cross  section  for  the  triplet  states  are 
discussed  in  the  next  section. 

The  rate  coefficients  for  all  the  Inelastic  processes  discussed  above 
have  been  obtained14*  15  for  a  Maxwellian  electron  velocity  distribution  and 
are  utilized  in  the  breakdown  code.  Thus,  if  one  designates  x^  as  the 
electron  impact  rate-coefficient  for  the  excitation  of  the  ith  state  whose 
excitation  energy  is  E^,  the  rate  of  the  electron  energy  loss  in  N2  can  be 
expressed  as 


Where  the  simulation  is  over  eight  ground  state  vibrational  levels,  the 
singlet  and  triplet  electronic  states,  the  dissociation  and  the  ionization 
of  Na-  A  similar  equation  for  the  energy  loss  in  the  nitrogen  atom  can  be 
utilized,  with  the  appropriate  rate  coefficients. 

Other  electron  energy  loss  processes  are  the  excitation  of  the  excited 
states  from  a  lower  state  (not  the  ground  state)  and  the  ionization  of  the 
excited  states. 

The  energy  gain  mechanisms  for  an  electron  are  the  energy  absorption 
from  the  field  of  the  microwave  i.e.,  the  inverse  bremsstrahlung,  the 
electron  impact  deexcitation  of  the  ground  state  vibrational  levels  and  the 
deexcitation  of  the  electronic  states  of  the  N2  molecule  and  the  nitrogen 

dtOID* 

3.3  EXCITED  STATES 

The  code  as  indicated  earlier  calculates  the  time  histories  of  several 
excited  states  of  Ns  which  are  discussed  below. 

A3£'-  The  A3£  state  is  a  metastable  and  has  a  life  time  of  1  sec  and  an 


excitation16  energy  of  6.1  eV.  Its  excitation  cross  section  by  electron 
impact  has  been  measured  and  calculated  by  several  investigators  which  are 
discussed  in  Ref  (7) »  and  shown  in  Figure  5*  The  electron  impact  ionization 
of  the  A3!  state  has  been  calculated  bu  Kukulin,  et  al.17  The  excitation 
energy  of  the  A3^  state  generally  ends  up  as  a  heating  source  for  the  gas 
because  of  the  rapid  quenching  of  the  state  by  the  atomic  nitrogen.  The 
rate  coefficient  for  the  quenching16  of  the  A3E  by  N  is  5  x  10-11  cm3/sec. 

B3!:  The  B3!  state  is  excited,  from  the  ground  state  of  Ng,  by  electron 
impact,  has  an  excitation  energy16  of  7-2  eV,  and  is  the  upper  level  for  the 
first  positive  band  system,  (B3tt  -♦  A3tt)  .  The  strong  vibrational  bands  of 
this  system18  are  given  in  Table  1. 


TABLE  1 


Relatively  Strong  Emissions  in  the  First  Positive  Band  System 


Wave  length  (A)  Transition  6j.  v) 


7503 

2 

6704 

5, 

2 

6623 

6, 

3 

6  544 

7, 

4 

6468 

a. 

5 

6394 

9, 

6 

6322 

10, 

7 

Wavelength  (l)  Transition  (6.  \j) 


6o69 

6, 

2 

6013 

7, 

3 

5959 

8, 

4 

5906 

9, 

5 

5854 

10, 

6 

5804 

11, 

7 

5755 

12, 

3 

The  life  times  of  these  transitions  are  in  the  range  of  10“  5  sec  to 
10"6  sec  and  are  quenched19  by  N2  with  a  rate  coefficient  of  ~  (1-2)  x  10  “11 
cm3/ sec. 

The  electron  imoact  excitation  cross  section  for  the  B3rr  state  have  been 
measured  and  calculated  by  numerous  investigators7  and  some  of  these  are 
shown  in  Figure  6. 

C3tt:  The  C3*r  state  is  excited  from  the  ground  state  of  Ng,  by  electron 
Impact  has  an  excitation  energy16  of  11  eV  and  is  the  upper  level  for  the 
second  positive  band  systems,  (C^tt  -♦  B3tt)  .  Some  of  the  strong18  bands  of 
this  system  are  given  in  Table  2. 
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TABLE  2 


Relatively  Strong  Emissions  in  the  Second  Positive  Band  System 
Wavelength  Transition  (•&.  0)  Wavelength  (l)  Transition  (i.  <!i) 


4059 

0,  3 

3755 

1.  3 

3998 

1,  ^ 

3710 

A 

OJ 

3943 

2,  5 

3577 

0,  1 

3894 

3,  6 

3536 

1,  2 

3804 

0,  2 

3371 

O 

A 

O 

3159 

1,  0 

3163 

2,  1 

The  lifetimes  of  these  transitions  are20  in  the  range  of  tens  of 
nanoseconds  to  fractions  of  a  microsecond.  Those  vibrational  levels  are 
quenched  by  N2  where  the  rate  coefficient  for  the  quenching  of  the  (0,  0) 
transition  is21  1.15  x  10” 11  cm3/ sec. 

The  electron  impact  excitation  cross  section  for  the  second  positive 
band  system  has  been  calculated  and  measured  extensively,7  especially  the 
cross  section  for  the  excitation  of  the  (0,  0)  transition  at  3371  A.  The 
excitation  cross  section  is  shown  in  Figure  7. 

N2+(B) :  The  nitrogen  molecule  has  several  ionization  contlnua  where 
upon  the  ionization  the  molecular  ion  is  in  an  excited  state.  These  are 
N2+(A),  Na+ (B) ,  Na+(C),  etc.  The  Na+(B)  state  has  an  ionization  threshold  of 
18. 8  eV  and  is  the  upper  level  for  the  first  negative  band  system, 

N2  (B2^  -♦  X5^;).  The  strong18  emissions  from  this  band  system  are  given  in 
Table  3. 
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TABLE  3 


Relatively  Strong  Emissions  in 

the  First  Negative 

Band  System 

Wavelenath  (k) 

Transition  (4.  G) 

Wavelength  (A) 

Transition  (u.  \i) 

5228 

0,  3 

4236 

(1,  2) 

5076 

2,  5 

3914 

(0,  0) 

4709 

0,  2 

4599 

2,  4 

4515 

4,  6 

4278 

0,  1 

The  lifetimes  of  these  states  range33  from  tens  of  nanoseconds  to 
microseconds  and  are  quenched  rapidly  by  N2.  The  rate  coefficient22  for  the 
quenching  of  the  zeroth  vibrational  level  is  4  x  10  _lD  cm3/ sec. 

The  fractional  ionization  leading  to  the  Na+ (B)  state  is  well 
established'  for  ionization  by  electrons  with  energy  of  100  eV  and  higher. 
However,  some  data  exist23  for  the  fractional  ionization  below  100  eV,  and 
emission  cross  section34  which  should  be  of  more  relevance  to  the  ionization 
of  N2  by  the  microwave  radiation. 

NC^):  The  metastable  state  of  the  nitrogen  atom  has  an  excitation 
energy16  of  2.37  eV  and  arises  by  the  dissociation  of  N2*  the  dissociative 
recombination  of  N2  ,  and  by  direct  electron  impact  excitation  from  the 
ground  state  of  N.  The  electron  impact  excitation  cross  section  of  N(®D) 
has  been  calculated.7 

3.4  GAS  HEATING 

One  of  the  basic  interests  in  the  air  breakdown  is  the  heating  of  the 
plasma  and  the  neutral  species.  The  heating  of  the  neutral  species  and  the 
ions  (since  the  collision  between  heavy  particles  is  very  efficient  we 
assume  that  the  ions  and  the  neutrals  have  the  same  temperature)  occurs 
through  several  atomic  and  molecular  processes.  These  are:  the  elastic 
collisions  of  the  electrons  with  the  heavy  particles,  the  quenching  of  the 
excited  electronic  states  by  the  neutral  species,  the  dissociative  recombina¬ 
tion  of  the  molecular  ions  and  the  charge  exchange  process.  The  quenching  of 
the  vibrational  energy  through  the  collision  between  neutral  species  and 
vibrationally  excited  molecules  is  another  essential  source  for  the  heating 
of  the  gas.  However,  the  slowest  heating  rate  among  these  various  processes 
is  the  vibrational  -  translational  energy  exchange. 


18 


When  Che  gas  temperature  rises,  two  additional  ionization  processes 
become  operative.  These  are  the  thermal  ionization  and  the  associative 
ionization  as  expressed  by  Equations  (20)  and  (21),  respectively. 

N2  +  N2  -»  N2+  +  N2  +  e  (20) 

N  +  N  -*  N2+  +  e  (21) 

4.  APPLICATIONS  TO  BREAKDOWN  STUDIES 

There  exist  considerable  experimental  25*26  data  on  the  microwave 
breakdown  in  air  and  other  gaseous  elements.  With  the  advent  of  lasers,  the 
air  breakdown  studies  at  optical  frequencies  have  increased  considerably?7-29 
Theoretical  calculations  for  the  air,  or  gas  breakdown  generally  proceed  from 
the  Boltzman's  equation  which  is  solved  self  consistently  for  the  electron 
energy  and  its  distribution.  Analytical  me thods2S» 27,30  are  also  utilized  to 
obtain  the  breakdown  threshold  power.  In  this  report  we  have  assumed  a 
maxwellian  electron  velocity  distribution  in  our  codes  (we  are  planning  to 
develop  a  Boltzman  code  for  this  purpose)  and  have  utilized  the  code  for  the 
breakdown  calculations. 

If.l  MICROWAVE  BREAKDOWN  IN  N2 

In  order  to  compare  our  results  with  experiments,  we  have  incorporated 

the  effect  of  diffusion  into  the  code.  Thus  we  define  an  effective  ionization 

rate  -  Da  to  allow  for  diffusion  in  our  code.  The  diffusion  is 

assumed  to  be  free  with  a  diffusion  coefficient  D  =“»  V2  t2  5  x  1015  T  . 

e  c  e 

The  time  t  *  1/v  is  the  electron-neutral  (electron-ion)  collision  time  and 
c  c 

Vg  is  the  electron  thermal  velocity.  The  diffusion  length  A  is  regarded  as 
being  specified  by  the  cavity  and  is  assumed  known.  Our  definition  of  break¬ 
down  is  defined  as  the  electric  field  required  to  achieve  a  net  increase  in 
the  electron  density  for  a  ten  nanosecond  pulse.  Once  the  cascade  process 
begins,  the  electron  density  quickly  reaches  the  critical  density.  The 
dominant  loss  processes  at  low  and  high  pressures  (Pal  Torr)  are  diffusion 
and  dissociative  recombination  (N4+  +  e  ->  N2  +  N2) ,  respectively.  Our 
comparison  will  be  to  experimentally  measure  cw  breakdown  in  nitrogen. 

The  criterion  for  cw  breakdown  in  nitrogen  is  M  vr  +  D/A2  where  v  is  the 

dissociative  recombination  rate31,  v  ■  5»0  x  10“S  N4+/T  ,  where  T  is  in  eV. 

r  e  e 

Both  the  diffusion  length  A  and  the  ambient  electron  density  (■%  -  10z  cm"3) 
were  specified  in  the  experiment  and  are  known  parameters  in  our  calculation. 


Calculations  were  performed  for  various  frequencies,  diffusion  lengths, 
and  pressures.  Figure  (8)  shows  the  results  obtained  by  MacDonald  et  al.26 
for  cw  breakdown  fields  in  nitrogen  at  frequencies  992MH  and  9.4GH  with 
diffusion  lengths  of  0.631  cm  and  0.4  cm,  respectively.  The  stars  in 
Figure  (8)  denote  the  calculated  breakdown  threshold  intensity  using  N2-Chem. 
The  calculated  thresholds  are  higher  as  expected  for  pulsed  breakdown.  Table 
4  shows  the  electron  temperature,  breakdown  intensity,  and  characteristic 
time  for  electron  production  t  -  (dN^/dt)  1  as  a  function  of  the  back¬ 
ground  pressure  at  f  *  922  MH  and  A  =  0.631  cm. 

Table  4 

Pulse  Breakdown  Threshold  Power  for  10  cm  Wavelength  for  Various  Nitrogen 

Pressures 


Neutral_ 
densitv  (cm  3) 

Electron 
Temnerature  (eV) 

Pulse  Breakdown 
Power  (Watts/ cm2') 

t  (sec) 

1.32  x  1018 

1-7 

2200. 

7.1  x  10‘7 

1.32  x  1017 

2.4 

80. 

4.3  x  10-7 

1.32  x  1016 

6.4 

o 

VO 

2.5  x  10‘7 

1.32  x  1015 

41.0 

153. 

3.7  x  10"8 

The  characteristic  time  r  represents  roughly  the  maximum  growth  time  for  the 
electrons  at  the  start  of  the  avalanche  process.  Experiments  show  that  once 
the  pulse  length  tp  j.  25ns,  the  pulse  breakdown  threshold  approaches  the  cw 
breakdown  limit.  The  generation  time  is  approximately  tg  *  j  log  2  and  we 
require  20  -  25  generations  to  reach  electron  densities  of  order 

~  1011  cm"3.  Again  our  definition  of  the  breakdown  power  is  the  intensity 
required  to  achieve  a  net  increase  in  the  electron  density  in  10  ns. 
Therefore,  our  results  are  certainly  applicable  to  pulse  breakdown  when 
tg  «  25p,3.  Since  this  is  easily  satisfied  by  our  results,  we  would  expect 
the  higher  breakdown  field  as  characteristic  of  pulse  breakdown.  The  high 
electron  temperature  at  low  pressure  results  from  the  strong  diffusive  losses, 
The  agreement  at  lower  pressures  may  be  somewhat  fortuitous  since  we 
extrapolate  the  ionization  and  excitation  rates  to  energies  above  20  eV. 
Current  rate  coefficients  in  the  code  are  tabulated  for  T^  up  to  20  eV. 


20 


i  I  n  null _ I  l  mill] _ 1  ]  ]  LI  Jill _ I  1  1 

0.1  1.0  IOO  too 

PRESSURE  (TORR) 


Fig.  8  —  Comparison  between  experimental  CW  N2  breakdown  threshold  power 
and  calculations  using  the  Nj  chemistry  code,  CHEMN2.  Calculations  are 
indicated  by  circles  and  crosses  for  pulsed  and  CW  breakdown  fields,  respec¬ 
tively. 


However,  data  is  available  for  T  above  this  temperature  and  will  be 

e 

incorporated  into  the  code  if  desired, 

k.2  MICROWAVE  AIR  BREAKDOWN  OFF  A  REFLECTING  SURFACE 

We  have  also  coupled  the  E  and  M  code  to  a  detailed  air  chemistry 
code32  (CHMAIR)  to  study  the  breakdown  in  air.  We  have,  however,  chosen  to 
calculate  the  breakdown  off  a  reflecting  surface  to  illustrate  the  breakdown 
in  air  as  a  result  of  the  interference  between  the  incident  and  the  reflected 
waves.  In  these  calculations  we  neglect  the  hydrodynamic  effects  and  assume 
that  chemical  processes  are  more  effective  in  depleting  the  electron  density 
than  diffusion.  This  is  especially  true  in  microwave  absorption  experiments 
with  short  pulses  (t  ^  40  ns)  and  rather  long  wavelength  radiation  conducted 
at  atmospheric  pressures.  Assuming  the  diffusion  length  A  ="  X»  ft 
atmospheric  condition  the  characteristic  electron  diffusion  time,  for  10  cm 
radiation,  is  roughly  t^  =■*  (A/L^p)2  t£  =*  10  7  sec  where  is  the  mean 

free  path  of  the  electron.  Since  the  diffusion  time  tQ  »  t  =“  kOns,  our 
neglect  of  hydrodynamic  effect  would  seem  a  good  assumption. 

The  intensity  of  the  standing  wave  patterns  resulting  from  the 
interference  maximum  can  easily  be  above  the  threshold  for  air  breakdown. 

We  imagine  the  microwave  to  be  incident  on  a  reflecting  surface  with 
intensity  Iq  <  Ig.  At  the  first  maximum  (one-quarter  wavelength  from  the 
surface)  the  field  intensity  is  I  =  4  I,  with  air  breakdown  commencing 
rapidly  a  distance  -p“  from  the  surface  (Figure  9a) .  Figure  9b  shows  the 
rapid  increase  of  the  electron  density  as  the  avalanche  begins.  The  extreme 
sensitivity  of  the  air  breakdown  to  the  electric  field  essentially  localizes 
the  electron  density  resulting  in  sharp  scale  lengths,  L  =  Ng  (dN  /dy)“  1  <  \. 
In  this  simulation,  we  have  a  pulse  of  10cm  wavelength  radiation  incident 
normal  to  a  reflecting  surface  with  an  intensity  of  Iq  =  0.5MW/cm2  and  a 
pulse  duration  of  kOns.  The  initial  ambient  electron  density  is  taken  to  be 
tv,  =  10cm”3  at  standard  temperature  and  pressure  (STP) . 

In  a  collisional  plasma,  »  oj,  the  cut  off  density  for  propagation  of 
electromagnetic  radiation  is  not  sharply  defined  as  in  collisionless  plasmas. 
At  sufficiently  high  density  such  that  the  skin  depth  5  s  X,  the  breakdown 
region  closest  to  the  surface  is  decoupled  from  the  microwave  radiation. 

This  is  seen  in  Figure  10a,  where  we  show  the  square  of  the  electric  field  at 
the  end  of  a  kOns  pulse.  In  this  particular  simulation,  we  allow  for 
multiple  sources  (each  source  below  threshold)  such  that  the  net  intensity 
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Fig.  9  —  (a)  Air  breakdown  off  a  reflecting  surface  at  atmospheric  conditions 
with  incident  power  I  ■  0.5MW/cm2  at  A  =  10  cm.  (b)  Time  development  of 
the  electron  density  (Ne/Nc). 


(a)  (b) 


Fig.  10  —  (a)  The  square  of  the  electric  field  at  the  end  of  the  pulse  for  tp 
*  40  nsec,  (b)  Time  development  of  the  electron  density  (Ne/Nc). 
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Is  *  2MW/cm3  >  I0.  The  corresponding  density  for  6  <  \  can  be  easily 
calculated  to  be  N^  =“  JOON^.  Figure  10b  shows  the  time  history  of  the  rise 
of  the  electron  density  to  Ng  <  500Nc>  Though  air  breakdown  is  occurring 
throughout  the  volume  for  this  particular  simulation,  the  extreme  sensitivity 
of  the  air  breakdown  to  the  electric  field  again  localizes  the  region  of  the 
electron  density.  The  important  point  to  be  made  here  is  the  high  electron 
density.  As  will  be  shown  later,  the  reflection  coefficient  for  an 
idealized  sharp  electron  density  profile  is  insignificant  until  Ng  =“  lO5^  at 
\  =  10cm.  The  high  electron  density  makes  real  the  possibility  of  driving  a 
hydrodynamic  response  with  microwaves,  pending  the  inclusion  of  the 
hydrodynamic  effects  and  various  pulse  shapes  and  pulse  lengths.  Early  work 
by  Wood33  using  the  hydrodynamic  code  LASNEX  suggests  the  formation  of  shocks 
at  atmospheric  conditions  using  microwave  radiation. 

4.3  REFLECTION  OF  MICROWAVE  RADIATION 

Once  the  electron  density  becomes  sufficiently  large  N^  »  N^,  the 
plasma  becomes  a  good  absorber  of  microwave  radiation.  However,  at  too  high 
densities,  the  plasma  becomes  a  good  reflector  of  microwave.  We  wish  to 
obtain  maximum  bounds  on  the  reflection  and  absorption  coefficient  of  the 
plasma  as  a  function  of  wavelength  and  electron  density.  We  proceed  further 
by  calculating  the  absorption  and  reflection  coefficient  for  an  idealized 
slab  density  profile  of  density  N  and  thickness  d  (See  Figure  11).  We 


consider  multiple  reflection  and  vary  the  wavelength,  density,  and  thickness 
of  the  microwave  plasma  system.  In  particular  we  examine  the  reflection  at 
two  microwave  frequency  of  3GH  and  30GH,  respectively.  The  electron-neutral 


collision  frequency  is  assumed  to  be  vc  3  5  x  1012  sec"1  which  is  typical  at 


ambient  atmospheric  conditions.  The  reflection  and  transmission  coefficients 


are  given  by 

R  ■ 

Eref 

m 

Einc 

5  (1  -  e’ie) 

~1-  S*  e-1C 


E_ 

Trans 

2 

(1  -  S2)  e'le/2 

E. 

inc 

1  -  Sa  e‘ie 
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SIMPLE  PHYSICAL  MODEL 
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Fig.  11  —  The  idealized  slab  density  profile  is  characteristic  of  systems 
with  high  electron  densities  and  steep  density  profiles. 
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where  S  *  (1  -  n)/(l  +  n)  and  e  =  2n«)d/c.  The  complex  refraction  index 


n  =  n  -  in,  is 
r  i 


(1  -  a) 
2 


1  +  f  1  + 


kb‘ 


(1  -  a); 


* 


! 


(1  -  a) 


i+ 


4bJ 


(1  -  a)* 


Ji 


when  1  -  a  >  0  and 


when  1  -  a  <  0,  where  a 


n  -  1  ^  ‘ 

1) 

\Jl  +  _ 

4b2 
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■  (a 

-  D2  ' 1\ 

1  <a' 

1) 

\J1+  _ 

1  2 
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The  absorption  is  determined  by  the  condition  R  +  T  +  A  =  1  with  the 

thickness  measured  in  units  of  the  skin  depth  6>  5k  «  —  .  Since 

ni 

the  results  for  d  *  5  and  d  »  5  are  somewhat  similar,  we  present  the  results 
only  for  d  »  6. 

Figure  12  shows  the  reflection  and  absorption  coefficients,  as  a 
function  of  plasma  density  at  \  ®  10cm  and  \  *  1cm  light,  respectively.  If 


we  assume  an  absolute  electron  density  of 
coefficient  is  roughly  60%  for  10cm  light  (Ng 
light,  the  reflection  coefficient  is  approximately  25%  at  N 


1015cm*3 

•  lO4^  ). 


,  we  see  the  reflection 
Whereas,  for  1cm 


10 


15__-3 

cm 


102Nc.  Our  calculations  also  indicate  that  to  maintain  a  plasma  density  of 

N  « 
e 


1015cm"3  requires  an  absorbed  power  of  I  ^ 


0.8MW/cm2.  Using  a  10cm 
2  1.6MW/ cm2  >  I_. 

D 

Obviously,  once  the  incident  power  is  above  the  breakdown  threshold,  the 


light  would  require  an  incident  power  level  of  I.  a  1.6MW/cm2  >  I_. 

me  o 


plasma  would  be  impossible  to  maintain.  However,  the  prospects  of 
maintaining  a  plasma  at  Ng  ■  1015cm"  3  using  a  1cm  light  appears  feasible 
since  an  incident  intensity  of  I^nc  *•  lMW/cm2  <  1^  would  only  be  required. 
Though  the  electron  density  Ng  =-  1015cm  3  is  clearly  arbitrary,  it  is  not 

unreasonable  to  expect  the  plasma  to  become  a  good  reflector  when 

2  2  2 
OUp  >  vc  +  U)  . 

5.  Conclusions 

We  have  presented  and  discussed  two  codes  developed  for  the  studies  of 
the  pulsed  microwave  interaction  with  Na.  We  have  also  utilized  our  E  and  M 
code  coupled  to  a  chemistry  code  to  study  the  breakdown  in  air.  Our 
breakdown  threshold  power  calculations  agree  favorably  with  experimental  data 
in  N2. 

Our  calculations  for  breakdown  near  a  surface  brings  forth  the  importance 
of  reflection  from  a  plasma  on  its  maintenance  and  heating.  However, 
additional  studies  in  this  subject  are  in  progress  to  delineate  this  point 
further,  including  the  hydrodynamic  effects. 
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